Toroidal Crystals 
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Crystalline assemblages of identical sub-units packed together and elastically bent in the form of 
a torus have been found in the past ten years in a variety of systems of surprisingly different nature, 
such as viral capsids, self-assembled monolayers and carbon nanomaterials. In this Letter we analyze 
the structural properties of toroidal crystals and we provide a unified description based on the elastic 
theory of defects in curved geometries. We find ground states characterized by the presence of 5- 
fold disclinations on the exterior of the torus and 7-fold disclinations in the interior. The number 
of excess disclinations is controlled primarily by the aspect ratio of the torus, suggesting a novel 
mechanism for creating toroidal templates with precisely controlled valency via functionalization of 
the defect sites. 



Toroidal micelles can be formed in the self-assembly of 
dumbbell-shaped amphiphilic molecules [ij. Molecular 
dumbbells dissolved in a selective solvent self-assemble 
in aggregate structures determined by their amphiphilic 
character. This process yields coexisting spherical and 
open-ended cylindrical micelles which evolve slowly over 
the course of a week to more stable toroidal micelles. 
Toroidal geometries also occur in microbiology in the vi- 
ral capsid of the coronavirus torovirus . The torovirus 
is an RNA viral package of maximal diameter between 
120 and 140 nm and is surrounded, as other coronaviri- 
dae, by a double wreath/ring of cladding proteins. 

Carbon nanotori form another fascinating and techno- 
logically promising class of toroidal crystals 0] with re- 
markable magnetic and electronic properties. The inter- 
play between the ballistic motion of the tt electrons and 
the geometry of the embedding torus leads to a colossal 
paramagnetic moment in metallic carbon tori f3|. The 
high magnetic susceptibility, together with large ring ra- 
dius (up to 0.03 ^m for single tubes of 1.4 nm in diameter 
[^), implies that the mobility of the electrons in the high- 
est occupied molecular orbital is several times larger than 
that of other ring-shaped molecules such as benzene [y] . 

A unified theoretical framework to describe the struc- 
ture of toroidal crystals is provided by the elastic theory 
of defects in a curved background 0, i, i This 
formalism has the advantage of far fewer degrees of free- 
dom than a direct treatment of the microscopic interac- 
tions and allows one to explore the origin of the emergent 
symmetry observed and expected in toroidal crystals as 
the result of the interplay between defects and geometry. 
Since defective regions are natural places for biological 
activity and chemical linking, a thorough understanding 
of the surface topology of crystalline assemblages could 
represent a significant step toward a first-principle de- 
sign of entire libraries of nano and mesoscale components 
with precisely determined valency that could serve as the 
building blocks for mesomolecules or bulk materials via 
self-assembly or controlled fabrication. 

The embedding of an equal number of pentagonal and 
heptagonal disclination defects in a hexagonal carbon 
network was first proposed by Dunlap in 1992 as a pos- 




FIG. 1: (Color online) Isolated defects and scar phases in 
the (r, V) plane. When the number of vertices V increases 
the range of the screening curvature becomes smaller than 
one lattice spacing and disclinations appear delocalized in the 
form of a 5 — 7 — 5 grain boundary mini-scar. 



sible way to incorporate positive and negative Gaus- 
sian curvature into the cylindrical geometry of nan- 
otubes In the Dunlap construction the curvature is 
achieved by the insertion of "knees" in conjunction with 
each pentagon-heptagon pair arising from the junction of 
tubular segments of different chirality. The latter is con- 
ventionally specified by two integers {N, M) which iden- 
tify the direction along which a planar triangular lattice 
is rolled up in the form of a tube [l2l |. A junction between 
an {N, 0) and an (iV, N) tube is obtained, for instance, by 
placing a 7— fold disclination along the internal equator 
of the torus and a 5— fold disclination along the external 
3 3- By repeating the 5 — 7 connection pe- 



equator 

riodically it is possible to construct an infinite number 
of toroidal lattices with an even number of disclinations 
pairs and the dihedral symmetry group Dnh (where 2n 
is the total number of 5 — 7 pairs). 
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Another class of crystalline toroids with dihedral an- 
tiprismatic symmetry Dnd was initially proposed by Itoh 
at al [iBl shortly after Dunlap. Aimed at producing a 
structure similar to Cgo fuUerene, Itoh's original con- 
struction implied ten disclination pairs and the symmetry 
group D^d- In contrast to Dunlap toroids, the disclina- 
tions in the antiprismatic torus are not perfectly aligned 
along the equator but rather staggered at some angular 
distance Sip from the equatorial plane. A general clas- 
sification scheme for Dnd symmetric tori can be found 
in Ref. [l^. In this paper we will refer to the lattices 
themselves, with 2n disclination pairs, by the symbols 
TPn and TAn, with respective symmetry groups Dnh 
and Dnd- 

Within the elastic theory of defects on curved surfaces 
0, 0, the original interacting particle problem is 

mapped to a system of interacting disclination defects in 
a continuum elastic curved background. Disclinations are 
characterized by their topological or disclination charge, 
Qi, representing the departure of a vertex from a prefect 
triangular lattice. Thus qt = 6 — Ci, where Ci is the coor- 
dination number of the ith vertex. A classic theorem of 
Euler requires the total disclination charge of any trian- 
gulation of a two-dimensional manifold M to be equal to 
6xm, where xm is the Euler characteristic of M. In the 
case of the torus xm = 0, and thus disclinations must ap- 
pear in pairs of opposite disclination charge (i.e. 5— fold 
and 7— fold vertices with q,; = 1 and —1 respectively) in 
order to ensure disclination charge neutrality. The total 
free energy of a toroidal crystal with N disclinations can 
be expressed as 



N 



(1) 



where Y is the 2D- Young modulus and T{x) is the so- 
lution of the following Poisson problem with periodic 
boundary conditions: 



Ar(a;) = Yp[x) 



(2) 



where A is the Laplace-Beltrami operator on the torus 
and p{x) is the total topological charge density 
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(3) 



k=l 



of N disclinations located at the sites Xk plus a screen- 
ing contribution due to the Gaussian curvature K(x) of 
the embedding manifold. The first term in Eq. ([1]) rep- 
resents the long-range elastic distortion due to defects 
and curvature. Its form resembles the potential energy 
of a system of electrical charges. In this analogy the 
Gaussian curvature K{x) plays the role of a non-uniform 
background charge distribution while defects appear as 
positively and negatively charged point-like particles. As 
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FIG. 2: (Color online) Elastic energy of a 5 — 7 disclination 
dipole constrained to lie on the same meridian, as a function 
of the angular separation. In the inset, illustration of a circu- 
lar torus of radii Ri > 7?2. Regions of positive and negative 
Gaussian curvature have been shaded in red and blue respec- 
tively. A standard parametrization of the torus is obtained 
by considering the angles V' and 4>- 



a result, disclinations arrange themselves so to approxi- 
mately match the Gaussian curvature. The second term 
in Eq. ^ is the defect core-energy representing the en- 
ergy required to create a single disclination defect. This 
quantity is related to the short-distance cut-off of the 
elastic theory and is proportional to the square of the 
topological charge times a constant Cc- Finally is an 
offset corresponding to the free energy of a flat defect-free 
monolayer. 

Using standard analysis the function r(a;) can be writ- 
ten in the form 



r(.) ^ J 
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^^1k^d{x,Xk) - T,{x) , 
k=l 



(4) 



where rs{x) represents the stress field due to the curva- 
ture of the torus and is given by: 
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where r = R1/R2 is the aspect ratio of the torus. The 
function Td{x, Xk) is the stress field at the point x arising 
from the clastic distortion due to a defect at Xk and is 
given by 



Vd{x,Xk) K 
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where ^ = Ktan -'^(w tan -0/2) is a conformal angle aris- 
ing from the mapping of the torus to the periodic plane. 
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FIG. 3; (Color online) Phase diagram for Tp configurations 
in the plane {r,ec/AY). For r £ [3.68, 10.12] and ec ~ the 
structure is given by a Tio configuration with symmetry group 



z = ^ -\- icf), and a, k and uj are dimensionless constants 
depending on the aspect ratio r: 
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Finally -Qi is a Jacobian theta function and reflects the 
double periodicity of the torus [l3| . 

To analyze the elastic energy Eq. ([T]) we first consider 
the simplest case of two opposite sign disclinations lying 
on the same meridian of a torus of area A — 47r^i?ii?2- 
The elastic free energy of this system is shown in Fig. 
[2] as a function of the angular separation between the 
two disclinations. The minimum is obtained for the pos- 
itive disclination along the external equator of the torus 
(where K(x) is maximally positive) and the negative 
disclination along the internal equator (where K(x) is 
maximally negative). The picture emerging from this 
simple case suggests a procedure to systematically con- 
struct optimal defect patterns for an arbitrary number of 
disclination pairs, by placing the same number of equally 
spaced -1-1 and —1 disclinations along the internal and 
external equators respectively. We name this configura- 
tion with the symbol Tp, where p stands for the total 
number of disclination dipoles. 

A comparison between the free energy of different Tp 
configurations as a function of aspect ratio and core en- 
ergy is summarized in the phase diagram of Fig. [31 We 
stress here that only Tp graphs with p even have an em- 
bedding on the torus corresponding to lattices of the TP| 
class. Nevertheless a comparison with p— odd configura- 
tions can provide additional information regarding the 
stability of p— even lattices. The defect core energy en- 
tering Eq. II]) is equal to 2pec- Although dependent 



on the microscopic details of the system, the constant 
tcl-AY ^ 10~^ for a crystal of roughly 10'^ atoms. In 
the range r e [3.68, 10.12] and ec ^ 0, the structure is 
dominated by the Tio phase corresponding to a double 
ring of -1-1 and —1 disclinations distributed along the ex- 
ternal and internal equators of the torus at the vertices 
of a regular decagon. The corresponding lattice has D^h 
symmetry group. 

That this structure might represent a stable configura- 
tion for polygonal carbon toroids has been conjectured by 
the authors of Ref. [l3|, based on the argument that the 
36° angle arising from the insertion of ten pentagonal- 
heptagonal pairs into the lattice would optimize the ge- 
ometry of a nanotorus consistently with the structure 
of the sp^ bonds of the carbon network (unlike the 30° 
angle of the 6— fold symmetric configuration originally 
proposed by Dunlap). In later molecular dynamics simu- 
lations, Han [3l found that a 5— fold symmetric lattice, 
such as the one obtained from a (9,0)/(5,5) junction, is in 
fact stable for toroids with aspect ratio less then r ^ 10. 
The stability, in this case, results from the strain energy 
per atom being smaller than the binding energy of carbon 
atoms. We have shown here, from first principles, that a 
5— fold symmetric lattice indeed constitutes a minimum 
of the elastic energy for a broad range of aspect ratios 
and defect core energies. 

For small aspect ratios the 5— fold symmetric configu- 
ration becomes unstable and is replaced by the Tg phase. 
This configuration, however, does not correspond to a 
possible triangulation of the torus. In this regime, we 
expect the minimal energy structure to consist of ten 
non-coplanar disclination pairs as in the antiprismatic 
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FIG. 4: (Color online) Phase diagram of a 5— fold symmet- 
ric lattice in the plane {r,Sip). For small Sip and r in the 
range [3.3, 7.5] the prismatic TPS configuration is energeti- 
cally favored. For r < 3.3, the system undergoes a structural 
transition to the antiprismatic phase TA5. 
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TA5 lattice. The latter can be analyzed by introducing 
a further degree of freedom, S'ip, representing the angu- 
lar displacement of defects from the equatorial plane. A 
comparison between the TP5 and TA5 configurations is 
shown in Fig. U for different values of 6ip. For small 
(5-0 and r G [3.3, 7.5] the prismatic TPS configuration is 
energetically favored. For r < 3.3, however, the lattice 
undergoes a structural transition to the TA5 phase. For 
r > 7.5 the prismatic symmetry of the TP5 configuration 
breaks down again. In this regime, however, the elastic 
energy of both configurations rapidly becomes higher be- 
cause of the lower curvature and defects disappear. 

In the regime of large particle numbers, the amount 
of curvature required to screen the stress field of an iso- 
lated disclination in units of lattice spacing becomes too 
large and disclinations are unstable to grain boundary 
"scars" consisting of a linear array of tightly bound 5 — 7 
pairs radiating from an unpaired disclination 0, . In 
a manifold with variable Gaussian curvature this process 
yields the coexistence of isolated disclinations (in regions 
of large curvature) and scars. In the case of the torus 
the Gaussian curvature in the interior is always larger 
in magnitude than that in the exterior for any Ri > i?2 
and thus we may expect a regime in which the negative 
internal curvature is still large enough to support the 
existence of isolated 7— fold disclinations, while on the 
exterior of the torus disclinations are delocalized in the 
form of positively charged grain boundary scars. 

To check this hypothesis we compare the energy of the 
TP5 lattice previously described with that of "scarred" 
configurations obtained by decorating the original toroid 
in such a way that each pentagonal disclination on the 
external equator is replaced by a 5 — 7 — 5 scar. The 
result of this comparison is summarized in the phase di- 
agram of Fig. [l]in terms of r and the number of vertices 
of the triangular lattice V (the corresponding hexagonal 
lattice has twice the number of vertices, i.e. Vhex = 2^). 

V can be derived from the angular separation of neigh- 
boring disclinations in the same scar by approximating 

V ~ A/ Ay, with Ay = -^a^ the area of a hexagonal 
Voronoi cell and a the lattice spacing. When the as- 
pect ratio is increased from 1 to 6.8 the range of the 
curvature screening becomes shorter and the number of 
atoms required to destroy the stability of the TP5 lattice 
decreases. For r > 6.8, however, the geodesic distance 
between the two equators of the torus becomes too small 
and the repulsion between like-sign defects takes over. 
Thus the trend is inverted. 

On the material science side, defects can be functional- 
ized to provide binding sites for ligands. The number of 
excess disclinations will then determine the valency of the 
surface, which itself can serve as a building block for new 
molecules and bulk materials [l9| . The universality of 
the predictions presented here means that one has precise 



control of the valency by tuning the aspect ratio of the 
torus. This could lead to a very efficient scheme for cre- 
ating a variety of basic surfaces with well-defined valency 
which subsequently self-assemble into novel molecules 
and bulk structures. 
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